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First-principles calculations of the electronic structure of members of the i?NiC2 series are pre- 
sented, and their Fermi surfaces investigated for nesting propensities which might be linked to the 
charge-density waves exhibited by certain members of the series {R = Sm, Gd and Nd). Calcu- 
lations of the generalized susceptibility, xo(q, c^;), show strong peaks at the same q- vector in both 
the real and imaginary parts for these compounds. Moreover, this peak occurs at a wavevector 
which is very close to that experimentally observed in SmNiC2. In contrast, for LaNiC2 (which is a 
superconductor below 2.7K) as well as for ferromagnetic SmNiC2, there is no such sharp peak. This 
could explain the absence of a charge- density wave transition in the former, and the destruction of 
the charge- density wave that has been observed to accompany the onset of ferromagnetic order in 
the latter. 

PACS numbers: 71.45.Lr,71.18.+y 



I. INTRODUCTION 

The idea that the Fermi surface (FS), through instabil- 
ities in the electronic structure (such as nesting features 
or van Hove saddle points), can drive low-dimensional 
systems into new ground states such as spin-density or 
charge-density waves (CDW), is well-known (for exam- 
ple, the transition metal dichalcogenides [l] and Cr [2]). 
More recent theoretical studies of the role of the FS in 
CDW formation [3l [4] have contributed greatly to our 
understanding of the phenomenon, and have clarified a 
number of misconceptions. Specifically, those authors 
emphasize that simply inspecting the FS for possible 
nesting features without actually calculating the real part 
of the electronic susceptibility does not help in predicting 
CDW instabilities [4 . They further urge caution about 
attributing such nesting as the only (or even main) driv- 
ing force for such instabilities. In this paper we calculate 
the electronic structure, FS and susceptibilities of several 
members of the i?NiC2 family which have been shown 
to exhibit a fascinating array of electronic instabilities 
Oin]. SmNiC2 has recently been found to host an inter- 
esting interplay between CDW and ferromagnetic order 
[7^. On cooling below a temperature of 148K, a resistivity 
anomaly and the appearance of satellite peaks in x-ray 
scattering indicate the formation of a CDW. The criti- 
cal phonon softening, inferred from the thermal diffuse 
scattering above 148K (and which also disappears at the 
ferromagnetic transition) occurs at two specific wavevec- 
tors, namely qi = (0.5,0.52,0) and q{r = (0.5,0.5,0.5) 
[7]. Of these two wave vectors, it is the incommensu- 
rate qi modulation which develops into a CDW. Below 
17.7K the satellite peaks suddenly disappear and there is 
a sharp decrease in the resistivity, and both phenomena 
are coincident with the appearance of ferromagnetic or- 
der |7j. Shimomura et at. tentatitively suggest that the 
diffuse scattering associated with the phonon softening, 
the formation of the CDW state and their disappear- 
ance at the (first order) ferromagnetic transition could 
be collectively understood from a knowledge of the FS of 



SmNiC2 [7 . In this paper we address this issue through 
ah initio calculations of the electronic structure of the 
i?NiC2 system. 



II. ELECTRONIC STRUCTURE 
CALCULATIONS 

The i?NiC2 family possess the orthorhombic CeNiC2 
structure (Space Group 38, ^mm^)[8 , which is shown in 
Fig. [l^. The electronic band structure of i?NiC2 was 
calculated for several R (La, Nd, Sm and Gd) using 
the scalar-relativistic linear muffin-tin orbital (LMTO) 
method within the atomic sphere approximation (ASA) 
including combined-correction terms ^. Here, the po- 
tential is described within the local density approxi- 
mation for non-magnetic calculations [TO] and the lo- 




r X M Y r z 

^^^^^^^^^^ 




[eV] 

FIG. 1: (color online) The (a) crystal structure, (b) electronic 
band structure, and (c) density of states of SmNiC2. 
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La 
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4.563 
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4.34 


Nd 


3.765 


4.544 
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2.64 


Sm 


3.703 


4.529 


6.098 
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3.84 


Gd 


3.649 


4.518 


6.077 


Il2l 


4.28 



TABLE L Lattice parameters used in the electronic structure 
calculations, accompanied by the references from which they 
were taken. The density of states at the Fermi level, N{Ef), 
obtained from the calculations, is also given. 
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FIG. 2: (color online) The FS of SmNiC2, (a) in 3D with 
the first Brillouin zone marked, and (b) several contours in 
planes perpendicular to c*. The arrow in (b) denotes the 
experimentally observed CDW wavevector q = (0.5,0.52,0). 
All units are in terms of 27r/(a, 6, c). 



cal spin density approximation for spin-polarized calcula- 
tions [11 . The lattice constants used in the calculations 
are shown in Table [H All calculations included a basis 
of 5, p, d and (for Ni and R) f states; self-consistency 
was achieved at 1099 k-points in the irreducible (1/8)*^ 
wedge of the BZ (corresponding to a mesh of 13 x 13 x 13 
in the full BZ). The localized rare-earth /-electrons were 
described using the open-core method, in which they are 
treated as partially-filled core states. 

The electronic band structure and density of states 
of SmNiC2 are shown in Fig. and Fig. [T]^, respec- 
tively. The FS, shown in Fig.[2j comprises a single sheet, 
originating from hybridized states of predominantly rare- 
earth / nickel (i-character. A closer look at the FS reveals 
that it indeed has strong potential for nesting. In Fig.|2]3, 
a contour of the FS at several slices through constant 
c* is shown, with the arrow indicating a nesting vector 
of (0.5,0.52,0). In order to further investigate the role 
of nesting, the non-interacting susceptibility was investi- 
gated. Neglecting matrix elements, the non-interacting 
susceptibility, Xo(q, ^) for wavevector q and frequency uj 
can be expressed as: 



xo(q,^) 



/(en,k) - /(en',k+q) 



nn' ,k 



i5' 



(1) 



FIG. 3: (color online) Xo(q) for SmNiC2. The arrow de- 
notes the peak in both real and imaginary parts of Xo(q) at 
q = (0.5, 0.56, 0). 



where /(e) is the Fermi function and the sum runs over 
bands n and through the Brillouin zone. The Bril- 
louin zone integration is performed using the tetrahedron 
method of Rath and Freeman [M]. As recently empha- 
sized by Johannes et a/., although the imaginary part 
of Xo(q) depends only on details of the electronic struc- 
ture near the Fermi energy, it is the real part of Xo(q) 
(which has contributions coming from a bandwidth-sized 
window of energies) that could indicate a CDW instabil- 
ity [4 . Both the real and imaginary parts of the non- 
interacting susceptibility were calculated for SmNiC2, 
and are shown in Fig. |3] It is noteworthy that there 
is a strong peak at the same q vector in both the real 
and imaginary parts of Xo(<l)- Closer inspection reveals 
that the peak is at a wavevector of (0.5, 0.56, 0), which is 
very close to the CDW vector qi identified by Shimomura 
et al. [7 . It should be borne in mind when comparing 
the experimental qi and theoretical peak in Xo(q) that 
our calculations have been performed at the experimen- 
tal lattice constant, rather than at the energy minimum 
of the LMTO calculation. Nevertheless, excellent quali- 
tative, and good quantitative agreements are found be- 
tween experiment and these band structure calculations. 
Approximations beyond the ASA (such as full-potential 
schemes) and a fully-relativistic treatment may also im- 
prove the precise location of the peak in Xo(q)7 and its 
agreement with the experimental qi. The Xo(q) was also 
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FIG. 4: (color online) Real and imaginary parts of Xo(q) 
along the [1,1,1] direction, showing the peak at qr = 
(0.5, 0.5, 0.5) associated with the thermal diffuse scattering 
observed by Ref. [7 . Owing to the crystal structure, there is 
a symmetry axis at ^[i,i,i] = (0.5,0.5,0.5). 



calculated along the [1, 1, 1] direction (shown in Fig. [4| 
and again a peak was found at a vector of (0.5,0.5,0.5), 
which matches the qr vector associated with the thermal 
diffuse scattering by Shimomura et al [7 . 

Calculations for both NdNiC2 and GdNiC2 reveal sim- 
ilar electronic structure and Fermi surfaces to the Sm 
compound, and the nesting properties already identi- 
fied for SmNiC2 are still present in a visual inspection 
of the FS. Furthermore, calculations of both the real 
and imaginary parts of Xo(q) show comparable struc- 
ture, the strong peak previously identified for the Sm 
compound developing in i?e[x(q)] for q = (0.5,0.55,0) 
and q = (0.5,0.57,0) for NdNiC2 and GdNiC2, respec- 
tively (see Fig. [5|. Both the Nd and Gd compounds 
exhibit anomalies in the resistivity that have been iden- 
tified with the emergence of a CDW [5 . Although the 
wavevectors that define these CDW states have not yet 
been investigated, these results suggest the CDW struc- 
ture is expected to be similar in nature to that already 
observed for SmNiC2, with a weak dependence on the 
particular rare-earth compound (presumably associated 
with the lanthanide contraction). The coincidence of the 
peak structure in both the real and imaginary parts of 
the susceptilibity for several different rare-earths, demon- 
strated in Fig. |5j underlines the importance of the FS 
(through nesting) in these systems. 

Conversely, LaNiC2 shows no equivalent kink in its re- 
sistivity [5 for temperatures down to 12K, raising the 
question of whether this compound shares the same FS 
topology, and more specifically, whether the same nest- 
ing feature is present. We have addressed this by cal- 
culating the electronic structure and Xo(q) foi" LaNiC2. 
The density of states at the Fermi level is shown in Ta- 
ble [l[ alongside comparative values for the other rare- 
earth compounds addressed in this study. The FS, shown 



1.25- 



o 

in 
o 

II 



o 




0.75 



0.65 



FIG. 5: (color online) The real (top) and imaginary (bottom) 
parts of Xo(q) for different rare-earth elements (Nd, Sm and 
Gd), shown along the q = (0.5, ?7, 0) direction as a function of 
77. The dashed vertical lines indicate the location of the peaks 
in the real part of the susceptibility for each rare-earth. 



in Fig. [gJd, is found to be topologically quite different, 
and the features of the FS that give rise to the nesting in 
SmNiC2 appear strongly distorted, as demonstrated by 
the constant c* contours shown in Fig.joj^. The /m[xo(q)] 
shows a broad series of ripples along the q = (0.5,77,0) 
direction (see Fig. [t]), rather than the strong peak ob- 
served for SmNiC2, compounding the evidence of poor 
nesting. In the real part of the susceptibility, also shown 
in Fig.[7| the peak becomes spread out over some ~ 0.46* 
of the Brillouin zone, which is unlikely to drive a CDW. 
Calculations for SmNiC2 in the LaNiC2 structure yield 
similar results, emphasizing that the transition of the FS 
from strong to weaker nesting is a consequence of the 
structure, rather than the particular rare-earth substitu- 
tion. 

Spin-polarized calculations of the electronic structure 
have also been performed for SmNiC2 in order to investi- 
gate whether a similar transition in the nesting qualities 
of the FS (as a consequence of the spin-splitting of the 
FS) can explain the disappearance of the CDW state that 
accompanies the onset of ferromagnetism. These calcu- 
lations have been performed at the same lattice constant 
used for non-magnetic SmNiC2, since the structural tran- 
sition that has been reported at the Curie temperature 
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FIG. 6: (color online) The electronic band structure (a) and 
FS of LaNiC2, (b) in 3D with the first Brillouin zone marked, 
and (c) several contours as for Fig. |2] 



[5] is small and does not break any additional symmetry. 
Although a total ordered moment of 0.32/iB is reported 
experimentally [13], the magnetic state of the localized 
Sm / electrons has not yet been clarified. Measurements 
of X-ray Magnetic Circular Dichroism on the i^^-edge of 
Ni in SmNiC2 at 5K have indicated that there is a fi- 
nite moment on the Ni site [15 , but a full analysis of 
the moment was not possible. In our calculations, a spin 
moment of 0.17/iB/Sm develops, which is about half that 
reported experimentally, and is predominantly located 
on the Ni site. Inspection of the ferromagnetic FS sug- 
gests that the different band filling in the majority sheet 
substantially weakens the nesting, but it still remains to 
some extent in the minority FS. Calculations of Xo(<l) 
reveal a strong suppression of the (0.5, 77, 0) peak in the 
imaginary part, whereas its real counterpart displays a 
broad ridge over a large range of 77, which would again 
be unlikely to drive a CDW instability. Additional cal- 
culations at a fixed spin moment that corresponds to the 
(total) experimental magnetic moment reveal this effect 
is even stronger, and the broadening of the ridge struc- 
ture in the real part of x(q) is spread over almost half 
the Brillouin zone in the (0.5,7^,0) direction. 



III. CONCLUSION 

We have calculated the FS of several (La, Nd, Sm and 
Gd) members of the i?NiC2 series, accompanied by cal- 
culations of both real and imaginary parts of the non- 
interacting susceptibility Xo(q)- For the Nd, Sm and Gd 
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FIG. 7: (color online) Real (top) and imaginary (bottom) 
parts of Xo(q) for LaNiC2 and ferromagnetic SmNiC2 (at the 
LMTO moment of 0.17/iB)- For comparison, the Xo(q) for 
non-magnetic SmNiC2 shown in Fig. |3] is reproduced here. 
Note the broad hump in Xo(q) for both LaNiC2 and ferro- 
magnetic SmNiC2, in contrast to the sharp peak shown in 
Fig. [3] for non-magnetic SmNiC2. 



members, peaks in the imaginary part of xo, that can be 
directly identified with corresponding nesting properties 
of the FS, still persist in the real part of xo, which is the 
relevant quantity for a CDW instability. The wavevector 
[0.5,7^,0] corresponding to these peaks in Xo(q) is close 
[j] = 0.56) to the experimental wavevector (77 = 0.52) ob- 
tained from x-ray scattering measurements of SmNiC2 
[7 . On the other hand, the FS of LaNiC2 is found 
to be markedly different, and the Xo(q) does not ex- 
hibit such strong structure, providing a possible route 
towards explaining the apparent absence of a CDW state 
for this particular compound. Spin-polarized calculations 
for SmNiC2 at two different magnetic moments (0.17 and 
0.32/iB) demonstrate the gradual destruction of the nest- 
ing properties of the FS with increasing moment. The 
strong peak in Xo(q) that is present at zero moment be- 
comes washed out as the spin-polarized bands separate, 
in agreement with the suggestion of Shimomura et al. 
that the ferromagnetic transition in this compound ad- 
versely affects the nesting properties of the FS and leads 
to a destruction of the CDW state jTj. 

In summary, our findings suggest that for the Sm, Gd 
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and Nd compounds, the FS does indeed nest at the CDW 
wavevector recently observed for SmNiC2, and more im- 
portantly, this strong divergence survives in the real part 
of the susceptibility. Furthermore, calculations of two 
isostructural compounds for which CDW formation has 
not been observed show substantially weaker divergences 
in their respective susceptibility, lending support to the 
idea that FS nesting plays an important role in deciding 
the fate of these systems. 

Note added during revision 

During the preparation of this manuscript, a recent 
study of the electronic structure of LaNiC2 has been re- 



ported [T6^, in which the authors have relaxed the inter- 
nal crystal structure. They go on to comment that the 
FS they obtain is in excellent agreement with our results. 
In response to their investigation, we have repeated the 
calculation of LaNiC2 for these parameters, and find no 
significant change in the FS. 
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